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ABSTRACT 

We present Near Infrared imaging polarimetry of three regions of massive star 
formation, G192.16 — 3.82, Cepheus A, and W42. In W42 we have discovered a 
new bipolar nebula located at the far side of the HII region behind the visible 
cluster of exciting stars. The axis of this new nebula is aligned with the magnetic 
field threading the entire cluster region. Polarization in the bipolar outflow neb- 
ulosity associated with G192.16 is consistent with a single illuminating source, 
too faint to be detected at 2 /im. Polarization in the reflection nebulosity asso- 
ciated with Ceph A requires more than one illuminating source, although HW2 
is clearly dominant. In all three objects, the magnetic field in the outflow at 
distances greater than ~ 0.2 pc is radial. In G192.16 the magnetic field geom- 
etry closer than ~ 0.2 pc to the embedded star appears chaotic. For G192.16 
the outflow is not aligned with the surrounding magnetic field, which lies in the 
galactic plane. In Ceph A, the outflow axis could be interpreted as being aligned 
with the galactic plane, but the magnetic field threading the region is not. Only 
in the case of W42 is the magnetic field threading the HII region aligned with 
the mean field in the surrounding galactic plane. 

Subject headings: ISM: magnetic fields — stars: winds, outflows — techniques: 
polarimetric 
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1. INTRODUCTION 

Molecular outflows are a common occurrence in the evolution of young stellar objects 
(YSOs) of all masses (Lada 1985; Bachiller 1996). The driving mechanism for these 
bipolar outflows is still a subject of considerable debate and most of our knowledge of this 
phenomenon is based on the study of low mass YSOs (Cabrit, Raga, & Gueth 1997; Shang, 
Shu, & Glassgold 1998). This is understandable because low mass YSOs are more common, 
closer on average, and often are less obscured by dust than deeply embedded, high mass 
YSOs. Currently, models invoking centrifugally driven winds are a popular mechanism for 
driving these outflows in low mass YSOs (e.g. Kudoh, Matsumoto, & Shibata 2003). In 
these models the magnetic field is twisted into a spiral pattern rising vertically off of the 
star's dense, rotating circumstellar disk. Material entrained in this field is forced up the 
vertical axis and produces the observed collimated outflow (Pudritz & Ouyed 1997; Kudoh 
& Shibata 1997). Konigl (1999) reviewed efforts to extend these magnetically driven wind 
models of bipolar outflows in low-mass stars to high mass YSO's. It is not clear if the 
massive outflows associated with very luminous YSOs are simply a scaled up version of 
these magnetically driven winds, or a completely different process. 

Determining the magnetic field geometry in massive bipolar outflows could provide 
important constraints on these models, but this is a difficult observational task. These 
outflows are most easily studied when the outflow axis lies in the plane of the sky and 
the star's equatorial disk is viewed nearly edge-on. This projection produces the classic 
bow-tie morphology seen in many sources and allows both poles of the outflow to be 
easily separated spatially on the sky. If dust grains in the bipolar lobes are aligned by the 
magnetic field in the outflow, then polarimetry of starlight shining through the outflow 
can measure the projected magnetic field geometry. If the pitch angle of a spiral magnetic 
field originating in the YSO accretion disk produces a very tight spiral, then the measured 
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field geometry will be perpendicular to the flow (Hodapp & Eiroa 1989). Polarimetry at 
infrared wavelengths was used by Jones & Amini (2003) and Itoh et al. (1999) to study the 
magnetic field geometry in the massive molecular outflow associated with DR21. Most of 
their observations were consistent with a magnetic field pointing radially outward, aligned 
with the axis of the outflow. 

A similar conclusion was reached by Hodapp (1990) using optical polarimetry of stars 
in the vicinity of Cepheus A. However, optical polarimetry can not penetrate deep into 
obscured regions, and is more likely to measure the larger scale magnetic field geometry in 
which the star formation region is embedded. Infrared (IR) polarimetry has the advantage 
of being able to penetrate the considerable dust extinction that typically accompanies 
deeply embedded massive YSOs (Hodapp & Eiroa 1989). This allows us to probe lines of 
sight more likely to pass through the denser portions of the outflow. Imaging polarimetry 
provides additional information on the reflection nebulosity produced by the bipolar 
outflow. Furthermore, at IR wavelengths, reflection nebulosity in heavily obscured regions 
close to the driving stellar source can be probed. 

Polarimetry at longer wavelengths where the aligned dust grains are in emission avoids 
the need for background field stars. However, mm- and sub-mm polarimetry of regions 
with one or more massive outflows are sensitive to the magnetic field geometry associated 
only with warm dust near the massive YSO's, not the more distant dust in the bipolar 
outflow. Greaves, Holland, & Ward-Thompson (1997) found a rough correlation between 
the position angle of the magnetic field in the warm dust centered on the YSO and the 
direction of the outflow. When the outflow lies in the plane of the sky, the magnetic field 
tends to be perpendicular to the outflow axis, just the opposite of what is found further out 
in the outflow (Jones & Amini 2003; Itoh et al. 1999). 

In this paper we employ the technique of imaging polarimetry at H (1.65 jim) and 
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K (2.2 jum) of three compact regions with massive star formation: G192.16 — 3.82, W42, 
and Cepheus A. G192.16 and Ceph A are well known outflow systems and Ceph A has 
had optical and near-IR polarimetry of surrounding stars (Hodapp & Eiroa 1989; Hodapp 
1990). W42 was not previously known to contain an outflow, but we have discovered a very 
reddened, elongated reflection nebula in the cluster indicative of a bipolar outflow. For all 
three systems we present IR polarimetry of stars within and nearby the central star forming 
region and imaging polarimetry of the reflection nebulosity associated with the outflow. 

2. OBSERVATIONS 

All observations were made with NSFCAM (Shure et al. 1994) in polarimetry mode 
on the NASA 1RTF 3-m telescope at a plate scale of 0'/3 per pixel, resulting in a field of 
view of 77" square. In polarimetry mode NSFCAM utilizes a rotating half-wave plate at the 
entrance window of the camera and a wire grid polarizer in the secondary filter wheel. This 
allows the observer to form images at polarization position angles of 0, 45, 90 and 135° on 
the sky. From these we compute Stokes quantities Q (Jo — -^90) and U (I45 — /135) for stars 
and nebulosity within the image. Details of the NSFCAM (+Polarimeter) technique are 
given in Jones (1997), Jones & Gehrz (2000), and Kelley et al. (2004). Table 1 summarizes 
observational details and parameters of our targets. 

Photometric and instrumental polarization calibrations were achieved by observing 
stars from the list of Elias et al. (1982). Unreddened early type stars in this list were 
assumed to have no polarization at the 0.1% level. Measurements of several of these 
stars indicate that the instrumental polarization is less than 0.1%. Because our technique 
requires a finite time between waveplate positions, sky and transmission fluctuations create 
systematic errors for single measurements at the ±0.2% level. These systematic effects are 
the primary source of error for bright objects. Position angle and polarization efficiency 
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calibrations were made by observing SI in p Oph, which was assumed to have P# = 3.9%, 
9 H = 28° and = 1.95%, 9 K = 28°. Additional checks on the position angle callibration 
were provided by observing BN in Orion and AFGL 2591. Photometric magnitudes were 
transfromed from the natural IRTF system (Tokunaga, Simons, & Vacca 2002) to the 
modern MKO system using the transformation relations described in Kelley et al. (2004). 

Our ability to measure the polarization for stars depends on three factors. If the star 
is very faint, it will be impossible to measure the fractional polarization, given the typical 
values of a few percent for interstellar polarization. If the star is very bright, but has a very 
low intrinsic polarization, it will also be impossible to measure the polarization. Finally, 
stars with images confused with nearby stars, stars in the negative (sky) beam, or stars 
embedded in bright nebulosity will also be impossible to accurately measure. In this paper 
we report polarimetry of stars only if the signal-to-noise for the fractional polarization 
(P/ap) is greater than 3. 

3. DISCUSSION 

3.1. Grain Alignment 

Interstellar polarization has long been attributed to the alignment of interstellar grains 
through some interaction with the magnetic field, such as a modified form of the classic 
Davis- Greenstein mechanism (Davis & Greenstein 1951). There is little doubt that in the 
diffuse ISM between stars magnetic alignment of some sort is the principle mechanism at 
work (for a review of grain alignment, see Jones (1996); Lazarian (2003)). This conclusion 
is based on the pattern of optical polarization vectors in the Milky Way (Heiles 2000) and 
comparisons between Faraday studies at radio frequencies and optical/IR polarimetry in 
external galaxies (Jones 1996). However, in specific locations not probed by these tests, 
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other mechanisms might be dominant. These include the original mechanical mechanism of 
Gold (1952) that used streaming motions of the gas as well as more advanced treatments 
(Lazarian 1995) and magneto-mechanical mechanisms such ambipolar diffusion (Roberge, 
Hanany, & Messinger 1995). 

Although there is no direct evidence for grain alignment by these mechanical 
mechanisms at this time, one location they might be expected to play a role is in the 
outflows from massive YSOs, which will certainly have significant streaming motions in the 
gas. The direction of grain alignment (long axis) for the original Gold mechanism is along 
the direction of the flow. This is perpendicular to the grain alignment we find based on 
polarimetry of background stars shining through the outflows discussed below. However, 
not all alternatives to the Davis- Greenstein mechanism align the grains in this manner. 
With the caveat that other alignment mechanisms may play a role in massive YSO outflows, 
in this paper we will assume magnetic grain alignment by some modified version of the 
Davis- Greenstein mechanism (Jones 1996), where the direction of the polarization vector 
for light in extinction corresponds to the magnetic field direction. 



3.2. G192.16-3.82 

G192.16 has been extensively studied in CO (Snell, Dickman, & Huang 1990; Shepherd, 
Watson, Sargent, & Churchwell 1998), other molecules (Molinari, Brand, Cesaroni, & 
Palla 1996), maser emission (Codella, Felli, & Natale 1996), and near-IR imaging (Hodapp 
1994). Molecular line maps show a single a bipolar outflow from a young star, most likely 
unconfused with other outflows, extending 1 — 2' to the East and West at a position angle 
close to 90°. There is an infrared reflection nebula that coincides with the molecular outflow 
(Hodapp 1994). The outflow is centered on IRAS 05553+1631, which is coincident with 
a cm- and mm- continuum source (Shepherd & Kurtz 1999; Shepherd, Claussen, & Kurtz 
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2002). The radio continuum fluxes are consistent with an early B star at a distance of 2 Kpc 
(Shepherd, Watson, Sargent, & Churchwell 1998). High angular resolution cm- and mm- 
continuum observations show evidence for a solar system size accretion disk surrounding the 
exiting star (Shepherd, Claussen, & Kurtz 2002). G192.16 has a very long bipolar outflow, 
about 3 pc, with Herbig-Haro objects extending up to 10 pc from the star (Devine et al. 
1999). 

Our K-band image of G192.16 core and portions of the outflow to the East and West 
is shown in Fig. 1. The .fT-band polarization vectors of six stars located along the line of 
sight are shown in the top panel of Fig. 1. The coordinates and polarization of these stars 
are listed in Table 2. The cm- and mm- continuum source believed to correspond to the 
central source of the outflow is located just to the West of the bright reflection nebulae in 
the central image, marked by a white cross. The white patches in Fig. 1 are due to stars in 
the negative (sky) beam. 

The polarization of stars seen toward G192.16 on the sky is shown in the larger context 
of the surrounding ISM in Fig. 2, where we include two stars (#7 and #8) well off the 
outflow axis to the South (also listed in Table 2). These two stars have polarization vectors 
closely aligned with the galactic plane at a position angle of 140°. For the distance and 
direction of G192.16, we are looking largely across spiral arms in the Milky Way and would 
expect the magnetic field to lie predominately in the plane (e.g. Jones 1996). Examination 
of this region of the sky in the large data base of optical polarimetry of stars compiled by 
Heiles (2000) also shows the interstellar polarization vectors lie in the plane of the galaxy. 
Thus, G192.16 is embedded in a large-scale magnetic field aligned with the galactic plane. 
The outflow axis of G192.16 is not aligned with the surrounding magnetic field, but is 
oriented at a significant angle of ~ 50° to the galactic plane. 

With one exception (#3), polarization vectors for the six stars in the same region of 
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the sky as the outflow are not aligned with the interstellar magnetic field. This strongly 
suggests the interstellar polarization of these stars is influenced by dust associated with 
G192.16. Stars #1 and #6 are 1' NE and SW of the core respectively and have 
polarization vectors pointing radially away from the central source. These stars are likely 
sampling the magnetic field geometry in the outflow, indicating that at a projected distance 
of 0.6 pc from the embedded YSO the magnetic field is largely radial, as was the case for 
DR21. However, within a projected distance of a few tenths of a parsec, the polarization 
vectors are more chaotic. Star #3 is aligned with the local ISM and may be a foreground 
object, but the other three stars show no simple pattern of polarization position angles. 
Unlike the DR21 outflow, where the observations were consistent with a purely radial field 
geometry everywhere in the outflow, the magnetic field geometry close to the exciting star 
in G192.16 must be more complex. 

Deep, high resolution imaging of the region surrounding the YSO at 2.1 /im has been 
reported by Indebetouw et al. (2003). They find a faint 2 fim source about 1" East of the 
radio and mm continuum source. The i^-band flux from this source shows little evidence 
for Br7 or shocked molecular hydrogen emission. The source is essentially undetected at 
H (1.65 /iin), indicating it is very red. Indebetouw et al. (2003) speculate that this source is 
either directly viewed photospheric emission from the central star or a bright knot of diffuse 
emission associated with the accretion disk. Below we examine the imaging polarimetry of 
this source and try and determine its true nature. 

fT-band polarimetry of the reflection nebulosity associated with the G192.16 outflow 
is illustrated in the bottom panel of Fig. 1 and listed in Table 3. These measurements 
were made using a synthetic 4.5" x 4.5" aperture with our /, Q and U images. All of the 
polarization vectors are consistent with simple scattering of light from a central illuminating 
source close to the location of the cm- and mm- continuum source. In Fig. 3 we examine 
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the region of bright nebulosity close to the illuminating source in more detail. The upper 
left-hand panel in Fig. 3 shows a grey scale image of the inner region. Here we see a bright 
patch of reflection nebulosity and a fainter source, indicated by the cross lines, about 3" 
West of the bright peak. This fainter source is the source discovered by Indebetouw et al. 
(2003) who had significantly better seeing (0.6") than in our image (1.7"). Based on our 
astrometry using POSS stars in the vicinity of G192.16, the faint source is located at RA = 
5 ft 58 m 13 s 6) Decl = + i6°3i'58" ± 1" J(2000), in excellent agreement with Indebetouw et al. 
(2003). 

Using a 1.5" x 1.5" synthetic aperture we made polarization measurements of the 
nebulosity shown in Fig. 3. The results of these measurements, uncorrected for seeing 
effects, are shown in the lower left panel of Fig. 3. The source discovered by Indebetouw et 
al. (2003) is highly polarized, and is very likely reflection nebulosity, not the central star 
itself. Because of seeing effects, flux from the brighter patch of nebulosity can contaminate 
the polarized intensity we measure for nearby regions. We can roughly subtract this 
contamination from patches nearby the bright source by using the point spread function 
of nearby stars to estimate the amount of flux underlying the patch of sky we are trying 
to measure. For the faint source, contamination from the bright patch was about 30% of 
the total intensity. When the correction is made to /, Q and U, we obtain the polarization 
vectors in the lower right panel of Fig. 3. These vectors point back to an illuminating 
source that is 1.5 ± 1.0" West of the suspected stellar source. Within our estimated position 
error of 1", this is also the location of the compact continuum source. We conclude that the 
massive YSO responsible for the molecular outflow and the illumination of the nebulosity 
in G192.16 is coincident with the cm- and mm- continuum source and is not visible at K. 

Our images are sufficiently deep to put new constraints on the extinction in front of 
the central star. In areas of our images well away from bright nebulosity, our 3 a detection 
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limit is about K = 18. Near the position of the exciting star, however, we are limited to 
about K = 16.5 due to the bright underlying nebulosity. We estimate the central star must 
be fainter than this limit. Using the expected value for the apparent K magnitude of a 
B2 star at 2 Kpc suggested by Indebetouw et al. (2003) of K = 9.7, we derive A K > 6.8. 
Using the extinction law from Rieke & Lebofsky (1985), this corresponds to greater than 
65 magnitudes of extinction at V. The infrared color for this amount of extinction is 
E(H — K) > 4.2, placing the H band brightness of this star well out of reach of our 
observations. 

3.3. W42 

W42 (= G25.4 — 0.2S) is a moderately compact, obscured HII region toward the inner 
Galaxy at I = 25.4°. Continuum radio maps at 5 GHz were made by Woodward, Heifer, 
& Pipher (1985) that showed an elongated emission structure with two cores smaller than 
the 6" beam. At the 5 Kpc distance they adopted, the radio continuum observations were 
consistent with an 06.5-7 star as the primary ionizing source. Blum, Conti, & Damineli 
(2000) have surveyed the stellar content of W42 with imaging and spectroscopy in the 
near-IR. They find a cluster of stars with a K band luminosity function similar to the 
Trapezium. The HII region is powered by a single star that they classify as 05-6.5, based 
on the classification system of Hanson, Conti, & Rieke (1997). Blum, Conti, & Damineli 
(2000) suggest a closer distance of ~ 3 kpc, based on spectroscopic parallax arguments 
applied to the O star and its dereddened K magnitude (see also Lester et al. 1985). The 
spectrum of other bright stars in the cluster are featureless, suggesting they are still YSO's. 

Our if-band image of W42 is shown in Fig. 4 along with the polarization vectors 
of fourteen stars located along a line of sight to the HII region. The coordinates and 
polarization data for these stars are listed in Table 4. Note that the polarization vectors 
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are relatively uniform in position angle. There is no evidence for more than one component 
to the magnetic field sampled along the line of sight to W42. The stars within 30" of the 
exciting O star (star #12 in Table 4, star #1 in Blum, Conti, & Damineli (2000)) have a 
mean position angle of 18° with a dispersion of 12°. The stars to the East and West of the 
central cluster have a position angle that lies in the plane of the Milky Way at about 28°. 
This is also the position angle of the long axis of the continuum radio emission observed by 
Woodward, Heifer, & Pipher (1985). 

Woodward, Heifer, & Pipher (1985) speculate that the elongated radio emission may 
be indicative of a bipolar source. Although there is no bipolar nebula clearly visible in our 
.fT-band image nor in the images in Blum, Conti, & Damineli (2000), close examination of 
the polarimetry does indeed reveal a bipolar nebula. In the left panel of Fig. 5 we show 
a grey scale image of the inner 30" of the HII region in total intensity in the K band. 
In the right panel we show an image in polarized intensity (Ip = \/ Q 2 + U 2 ) with the 
polarization vectors for five regions of highly polarized nebulosity superimposed. These 
five regions were measured using a synthetic aperture 2.7" x 2.7" in size. Highly polarized 
reflection nebulosity, almost unrecognizable in the total intensity image, is clearly visible in 
the polarized intensity image. The polarization vectors are consistent with an illuminating 
source located between the two brightest blobs of nebulosity, indicated by the cross marks 
on the right and bottom of the polarized intensity image. If there is a luminous YSO 
located between the two bright blobs of reflection nebulosity, it would have coordinates RA 
= 18 h 38 m 14!6, Decl. = -6°48'00" (J2000). The long axis of the reflection nebulosity is 
oriented at about 15°, essentially the same orientation as the magnetic field through the 
entire cluster. 

The newly discovered bipolar outflow is very heavily reddened by extinction. The 
H — K color of the two brightest blobs is H — K ~ +2.4, redder than any of the surrounding 
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stars (Fig. 4 in Blum, Conti, & Damineli 2000). Since reflection nebulosity can be 
intrinsically much bluer than H — K = 0.0 (depending on the spectral energy distribution 
of the illuminating star), our measured color translates to a visual extinction of Ay > 55. 
The illuminating star is presumably behind additional dust associated with an equatorial 
accretion disk, explaining why it is not visible at K. Blum, Conti, & Damineli (2000) argue 
that we see the association and its primary O star on the front side of a dense, molecular 
cloud. The newly discovered bipolar nebula is clearly deeper into this denser region than 
the cluster, indicating ongoing star formation just behind the interface between the HII 
region and the molecular cloud. This reminds us of the Trapezium, which has a comparable 
luminosity and a similar geometry with respect to our line of sight from the Earth. 

Optical polarimetry of stars in the surrounding sky taken from the compilation of 
Heiles (2000) clearly shows a general magnetic field direction in the diffuse ISM surrounding 
W42 that lies in the plane of the Galaxy at PA = 28°. The magnetic field threading W42 
is close to this position angle. We are presented with a picture of a very uniform magnetic 
field geometry threading through the entire cluster, and linking up smoothly with the field 
in the surrounding diffuse ISM. The newly discovered bipolar nebula has a long axis aligned 
in the same direction as the threaded magnetic field, indicating that the polar axis of the 
illuminating star is also aligned with the local magnetic field. This is in contrast to G192.16 
and DR21, where the axis of the bipolar nebula and the magnetic field in the outflow is well 
away from alignment with the galactic plane. 

3.4. Ceph A 

Cepheus A East contains a well studied, relatively nearby molecular outflow driven by 
one or more massive, very young stars. At a distance of 725 pc (Blaauw, Hiltner, & Johnson 
1959) the Far-IR luminosity is about 2.4 x 10 4 L Q , which is consistent with several early 
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B stars. The bipolar outflow is quite complex (Torrelles et al. 1993; Narayanan & Walker 
1996) and is suggestive of multiple outflow sources. The molecular line observations have 
also been interpreted as clumps in two sides of a poorly collimated bubble powered by one 
or more YSOs (Corcoran, Ray, & Mundt 1993). 

Radio continuum observations at a resolution of 1" (Hughes & Wouterloot 1984) reveal 
several compact sources, which they associate with early B stars. Source number 2 in their 
map, commonly referred to as HW2 in the subsequent literature, is usually considered to 
be the primary source of the outflow. This source was undetected at K by Casement & 
McLean (1996), who present imaging polarimetric observations of Ceph A on a 2.7" grid 
at J, H and K. Analysis of the highly polarized reflection nebulosity surrounding HW2 
allowed Casement & McLean (1996) to show that the origin of the scattered light from 
dust in the outflow is primarily from HW2. They also found that within 10 — 15" of HW2, 
the polarization vector pattern was more complex than simple single scattering and clearly 
required multiple scattering. 

Goetz et al. (1998) made a detailed study of Ceph A East using broad and narrow band 
imaging from 1.4 — 5 /im. They find two regions of shock-excited line emission, one 20" to 
the NE of HW2 and one 20" to the ESE of HW2. They identify HW3 (not HW2) as the 
primary driving source of the shock emission to the ESE and consequently favor a multiple 
outflow model. They find that the reflection nebulosity to the NE of HW2 and HW3 is 
exceptionally red in their K — L' image, indicating it is behind significant extinction. Goetz 
et al. (1998) were unable to detect HW2 at L' or M, implying A v > 200 if it is a bare Bl 
star. In their schematic diagram of the Ceph A region, Goetz et al. (1998) associate the 
rotating core observed in CS line emission by Narayanan & Walker (1996) with a 'Fluffy 
Torus' centered on HW2 and oriented at a position angle of about 135°. 

Glenn, Walker, & Young (1999) measured the polarization toward HW2 at 1.3 mm and 
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found the magnetic field was oriented at 170°, which is roughly perpendicular to the outflow 
axis. This result is consistent with the conclusion of Greaves, Holland, & Ward-Thompson 
(1997), who argue that when the bipolar outflow lies in the plane of the sky, the magnetic 
field observed at the Far-IR/sub-mm peak is perpendicular to the outflow axis. Hodapp 
(1990) observed I band polarimetry of nearby stars and found the magnetic field geometry 
V — 4' away from HW2 was inconsistent with a single component. Unlike G192.16 and W42, 
the optical polarimetry of stars within a few degrees of Ceph A does not present a clear 
picture of the surrounding interstellar magnetic field. Using the data base of interstellar 
polarimetry compiled by Heiles (2000), we plot the polarization vectors for stars within a 
few degrees of Ceph A in Fig. 7. South of the star forming region the polarization vectors 
are nicely aligned with the galactic plane, indicating the interstellar magnetic field lies in 
the plane along that line of sight, as expected. However, 1° North and West, some stars 
show interstellar polarization vectors that lie in the plane of the Galaxy, while others lie 
at a position angle of about 125°. We will argue below that the foreground interstellar 
magnetic field is aligned with the plane of the Milky Way disk (PA = 65°), the magnetic 
field threading the Ceph A complex is almost perpendicular to the plane (PA = 125°) and 
that the magnetic field in the outflow is radial with respect to HW2/3. 

A false color image of Ceph A using our H and K band images is shown in Fig. 6. The 
location of HW2 is shown as a green cross. The reddest tones correspond to H — K = +4.0 
and the bluest tones correspond to H — K — +1.75, which is still quite red. Several stars 
are bluer than H — K = +1.75 and they are saturated at a blue color. Depending on the 
intrinsic colors of the illuminating source, the H — K color of reflection nebulosity would 
usually be significantly bluer than any of the values we measure in Ceph A. The entire 
reflection nebulosity must be behind significant interstellar extinction, very likely associated 
with the star forming region itself. This intervening extinction is optically thickest close to 
the central source, and becomes thinner as one moves North and East from the massive 
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YSO. A similar conclusion was reached by Goetz et al. (1998) from their K — V image. 

In Fig. 8 we plot a number of features seen in Cepheus A over an underlaying grey scale 
of the K band intensity. The solid contours are from the radio continuum map of Hughes 
& Wouterloot (1984). Their source 2 (HW2) is associated with the primary central massive 
YSO that powers and illuminates the outflow. The narrow rectangles are the polarization 
vectors of the reflection nebulosity on a 2" grid. The solid black lines are the if -band 
polarization vectors of the five stars listed in Table 6. The position angle of the galactic 
plane is 65°, very close to the position angle of star #4. 

Based on the H — K colors we observe for the nebulosity, the extinction within 10 — 15" 
of HW2 is too great for background stars to shine through (see Goetz et al. 1998). Star $T 
in Table 6 is outside the bright nebulosity and is relatively red. It has a position angle very 
close to the position angles measured by Hodapp & Eiroa (1989) in H band for stars a few 
arc minutes north and northwest of Star #1. We hypothesize this star is either embedded 
within the Ceph A complex or less likely, shining through from behind. This would imply 
that the interstellar magnetic field in the immediate vicinity of Ceph A is at a position 
angle of ~ 125° (Hodapp & Eiroa 1989). 

Star #2 is much bluer than Star #1 and has a similar position angle, but is still 
outside the bright nebulosity and may also be sampling aligned grains in Ceph A. Star #3 
is much bluer than the underlying reflection nebulosity and can not be shining through the 
outflow. Polarimetry of Star #3 is sampling the magnetic field geometry well in front of the 
outflow and is consistent with a foreground field lying in the plane of the Galaxy. Stars #4 
and #5 are only slightly bluer than the underlying nebulosity and are probably sampling 
aligned grains associated with the near side of the Ceph A outflow. Polarimetry of Star 
#4 could be interpreted in two ways. It could be dominated by the foreground extinction, 
causing the position angle to lie in the plane. But, the position angle points directly back 
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to HW2/3 and could be indicative of a magnetic field in the outflow pointing radially to 
the exciting B stars. The position angle of star 5 also points directly back to HW2/3 and 
is most likely sampling the magnetic field in the outflow as well. Note that Hodapp (1990) 
and Hodapp & Eiroa (1989) measured several stars 1 — 2' East of Star #5, all of which have 
position angles very similar to Star #5. We argue that the magnetic field in the outflow at 
a distance greater than 20" from the cluster of massive YSOs, is radial. 

A polarization map at a spatial resolution of 0.9" is shown in Fig 9. The polarization 
vectors for the reflection nebulosity are dominated by single scattering directly from the 
central sources at distances of more than 5 — 10" from HW2/3. The polarization vectors to 
the North and West of the embedded stars are consistent with scattering of light only from 
HW2. The polarization vectors to the NE and E, however, are illuminated from a location 
about 5" SE of HW2. We do not see evidence for this at a coarser grid spacing in the 
polarization map of Casement & McLean (1996). Our polarization map requires more than 
one source of light being scattered by dust in the Ceph A outflow. Dust in the outflow to 
the N and W is illuminated only by the YSO associated with HW2. The dust in the outflow 
to the NE and E is seeing at least one additional source to the SE of HW2, perhaps source 
HW3. Goetz et al. (1998) associate the shocked gas coincident with HW7 in the radio map 
with outflow from HW3, not HW2. Note that both the axis of the elongated radio emission 
in HW7 and the polarization vector for star 5 point directly back to HW3. It may be that 
HW3 has brightened at near-IR wavelengths since the observations of Casement & McLean 
(1996) and is now contributing more to the reflection nebulosity. 

The polarization vectors strongly depart from simple centro-symmetric scattering 
along a band stretching across HW2 from the SE to the NW at a position angle of ~ 135°. 
This morphology corresponds well in position angle to the smaller 'Fluffy Torus' described 
in Goetz et al. (1998) and the rotating gas described by Narayanan & Walker (1996). 
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Casement & McLean (1996) conclude that this departure from centro-symmetric scattering 
as likely due to multiple scattering in dust near the central source. We are probably looking 
through an optically thick disk of material extending at least 15" (0.05 pc) to the NW and 
SE of the central cluster of B stars. Light passing through this region will be scattered at 
least twice and loose the simple centro-symmetric polarization pattern produced by single 
scattering. 

Based on all of the available observations of interstellar polarization toward Cepheus 
A, we propose the schematic picture of the magnetic field geometry within and in front of 
Ceph A shown in Figure 10. The foreground (and possibly background) magnetic field is 
closely aligned with the galactic plane at a PA of ~ 65°. The field threading the overall star 
forming complex is at a very different orientation of ~ 125°. At the Far-IR/sub-mm peak, 
essentially at the location of HW2 and HW3, the magnetic field is roughly perpendicular to 
the outflow at ~ 170°. We say roughly, because there are probably multiple outflows and 
the reflection nebulosity presents a very broad fan-like structure on the sky. Note, however, 
that the field direction at the Far-IR/sub-mm peak is perpendicular to the outflow axis in 
the model proposed by Narayanan & Walker (1996). Finally, the interstellar polarization 
to the East of the core cluster, as measured by Hodapp and ourselves, shows the magnetic 
field well out in the outflow is radial. 



4. CONCLUSIONS 

We have made imaging polarization observations of three regions of massive star 
formation at near-IR wavelengths. These observations allowed us to measure the interstellar 
polarization of stars along the line of sight to these regions and to measure the polarization 
of reflection nebulosity associated with embedded YSOs. Both G192 and Ceph A are 
well known bipolar molecular outflow sources with well identified obscured massive YSOs 
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driving the flows and illuminating the reflection nebulosity. In W42 we have discovered a 
previously unknown bipolar nebula and identified a probable location for the central star 
responsible for the observed reflection nebulosity. 

4.1. Reflection nebulosity 

The reflection nebulosity in G 192. 16 was bright enough in several locations for us 
to measure the polarization due to scattering of light from the central YSO by dust in 
the outflow. The observed polarization is entirely consistent with illumination by a single 
stellar source. This position of this object is coincident with the well studied point-like 
cm- and mm- continuum source previously identified as the illuminating YSO. The YSO is 
undetected in our maps at 2 jum down to K = +16.5 (3a). 

In Ceph A, the polarization map of the reflection nebulosity shows evidence for more 
than one source of illumination in addition to HW2, possibly HW3. Neither of these sources 
are detected at 2 /im in our maps. The polarization in the reflection nebulosity close to 
HW2 and stretching from HW2 for several arcseconds perpendicular to the outflow shows 
clear evidence of multiple scattering. This region probably contains the outer extremes of a 
large accretion disk surrounding HW2. 

4.2. The magnetic field geometry 

In W42 and Ceph A, the main axis of the outflow and the local magnetic field are 
within 15° of the galactic plane. In G192.16, the main axis of the outflow and the magnetic 
field in the outflow are misaligned by at least 60° with respect to the galactic plane. In 
W42 the magnetic field threading the HII region is well aligned with the galactic plane, and 
shows no evidence for additional components to the magnetic field geometry. The magnetic 
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field geometry within Ceph A is complicated, but we believe the field threading the complex 
is not aligned with the galactic plane, even though the outflow axis may be. In G192.16, 
the magnetic field geometry threading the compact HII region could not be determined. 

In all three regions, the local magnetic field is aligned with the direction of the outflow 
at distances of several tenths of a parsec and further from the central source. In G192.16 the 
magnetic field geometry closer than a few tenths of a parsec to the central source is probably 
chaotic, and does not exhibit a radial geometry. In Ceph A, we have no observations of 
the magnetic field geometry within 20" of the exciting star(s), but sub- mm observations 
indicate the field at the location of the embedded YSO(s) is perpendicular to the outflow 
axis. In W42, the newly discovered bipolar nebula is too small in angular extent on the 
sky for our technique to be able explore the embedded magnetic field geometry. The axis 
of the nebula is, however, exactly aligned with the magnetic field in the HII region and the 
galactic plane. 

The magnetically driven model for YSO outflows predicts a spiral magnetic field in 
the core of the outflow along the outflow axis. Our observations do not rule this model 
out if the field is a tight spiral only at small distances of less than ~ 0.2 pc from the 
star. At greater distances, our observations are consistent with a purely radial field. 
The tendency for the magnetic field at the location of the embedded YSO(s) (measured 
by sub-mm polarimetry) to lie perpendicular to the flow is probably due to a toroidal 
field in the extended equatorial disk, not in the outflow itself (see the various model 
geometries discussed in Greaves, Holland, & Ward-Thompson 1997). Neither our technique 
of measuring interstellar polarization of background stars nor sub-mm polarimetry with 
large beam sizes will be able to probe the magnetic field geometry in the outflow close to, 
but separated from, the embedded star and its accretion disk. High spatial resolution far-IR 
and mm- polarimetry will be necessary. 
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FOR HIGH QUALITY FIGURES VISIT http://www.astro.umn.edu/- tjj/ 
CLICK ON Outflow Paper Figures 

Fig. I. — Grey scale image of G192.16 — 3.84 at K (2.2 /im). The location of the mm- 
point source is indicated by a white cross. White objects in the frames are due to stars 
in the sky beam. The top panel has the polarization vectors for six of the stars listed in 
Table 2 superimposed on the image. The bottome panel has the polarization vectors for the 
reflection nebulosity measured at different locations in the outflow and listed in Table 3. All 
of the polarization vectors are consistent with light scattered from a single source located at 
the position of the mm- continuum source. 

Fig. 2. — Polarization of stars nearby and south of G192.16 — 3.84. The cross mark indicates 
the location of the mm- continuum source identified as the YSO powering the outflow. The 
molecular outflow is aligned roughly East- West. The galactic plane lays at a position angle 
of 120°. Note that the position angles of the two stars to the South (Table 2) show that 
the surrounding interstellar magnetic field lays in the galactic plane. However, the stars 
close to G192.16 are clearly at very different position angles than the magnetic field in the 
surrounding diffuse ISM. 

Fig. 3. — Detailed view of the reflection nebulosity near the central source in G192.16 — 3.84. 
The upper left panel shows a grey scale image at 2.2 /xm of the nebulosity. The cross marks 
indicate the location of the faint source discovered by Indebetouw et al. (2003). The thick 
solid line locates a cut through the image plotted in the upper right hand panel. This 
illustrates how flux from the bright peak will contaminate the fainter sources nearby due to 
seeing effects. The bottom left hand panel shows the polarization vectors at various locations 
in the region uncorrected for contamination by flux from the bright peak. The lower right 
hand panel shows these same polarization vectors after a rough correction for contamination 
from the bright peak. These vectors are consistent with a source of illumination 1.5 ± 1.0" 
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west of the source discovered by Indebetouw et al. (2003). 

Fig. 4. — Grey scale image of W42 (= G25.4 - 0.25) at K (2.2 fxm). Superimposed on the 
image are the polarization vectors for the stars listed in Table 4. The bright star in the 
center of the cluster (star 12 in Table 4) was identified by Blum, Conti, & Damineli (2000) 
as an 05-06.5 star and the primary ionizing source in the HII region. The surrounding 
interstellar magnetic field is closely aligned with the galactic plane. 

Fig. 5. — The left panel shows a grey scale image in total intensity at K of the inner 30" 
of W42. The bright star at center left is the O star identified by Blum, Conti, & Damineli 
(2000). The right panel shows the same region, but in polarized intensity. Note the highly 
polarized reflection nebulosity seen in the polarized intensity image is barely visible in the 
total intensity image. The polarization vectors for six locations in this reflection nebulosity 
are superimposed and the coordinates and polarimetry are listed in Table 5. The two arrows 
in the right hand panel indicate the location of a possible illuminating star for this reflection 
nebulosity. 

Fig. 6. — False color image of the Cepheus A outflow from H and K band images. The 
green cross indicates the location of the primary illuminating source. The reddest tones 
correspond to H — K = +4.0. The bluest tones correspond to H — K = 1.75. Stars with 
H — K colors bluer than +1.75 are saturated at the bluest color. The reflection nebulosity 
is clearly much more heavily reddened near the central star. 

Fig. 7. — Optical polarimetry of stars within a few degrees of Cepheus A (marked with a 
cross) taken from the compilation of Heiles (2000). The galactic plane is indicated by the 
dashed line. 

Fig. 8. — Grey scale image of Ceph A at K. The solid contours are the radio continuum 
observations of Hughes & Wouterloot (1984). Their Source 2 (HW2) is associated with the 
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central star powering and illuminating the outflow. The thin rectangles are the polarization 
vectors of the nebulosity measured on a 2" grid spacing. The polarization vectors of the five 
stars in Table 6 are plotted as solid black lines. 

Fig. 9. — Polarization vectors of the reflection nebulosity in Ceph A plotted on a 0.9" grid. 
The underlying intensity contours at K are shown to aid in location. The coordinates of 
HW2 are shown by a cross. 

Fig. 10. — Schematic diagram of the magnetic field geometry in Cepheus A, based on the 
interstellar polarization of stars and sub-mm polarimetry centered on the central cluster of 
massive YSOs. 



-28- 



Table 1. Observing Log 



Object 


Date 


Filters 


Mode 


G192.16 Core, East 


Feb. 14, 2002 


K 


Polarimetry 


G192.16 West 


Feb. 15, 2002 


K 


Polarimetry 


W42 


July 10, 2003 


H, K 


Polarimetry 


Cepheus A 


July 9, 2003 


K 


Polarimetry 


Cepheus A 


July 10, 2003 


H 


Polarimetry 



Table 2. Stars in G192.16 -3.84 



Star RAJ2000 Decl. P K (%) e P (%) 6 K {°) 

5 h 58 m +1 go 



1 


1K3 


32'30" 


1.6 


0.3 


52 


2 


1K6 


3219" 


2.3 


0.3 


93 


3 


14!6 


31'55" 


1.5 


0.3 


155 


4 


1K2 


31'39" 


1.3 


0.3 


87 


5 


1K7 


31'37" 


1.9 


0.3 


35 


6 


11!2 


31'21" 


2.6 


0.8 


178 


7 


14!3 


SO'OO" 


1.8 


0.4 


141 


8 


10:0 


30'04" 


2.2 


0.5 


138 
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Table 3. Reflection Nebulosity in G192.16 -3.84 



Patch RA J2000 Decl. P K (%) K (°) 
5 h 58 m +16° 



1 


18% 


31'51" 


45 


12 


2 




31'49" 


35 


12 


3 


1K4 


31'51" 


39 


15 


4 


13?8 


31'57" 


37 


11 


5 


12!9 


31'53" 


18 


165 


6 


1Z3 


32' 15" 


31 


46 


7 


1Z3 


32'00" 


17 


7 


8 


11?8 


31'39" 


62 


138 
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Table 4. Stars in W42 



Star RAJ2000 Decl. P K {%) e P (%) 6 K {°) 



1 


17:6 


47'21" 


0.9 


0.3 


32 


2 


16*3 


47' 15" 


2.2 


0.3 


5 


3 


17:7 


47'54" 


5.1 


0.9 


41 


4 


1&.2 


47'45" 


2.2 


0.3 


24 


5 


15!8 


47'45" 


2.1 


0.3 


178 


6 


15 S A 


47'39" 


0.8 


0.8 


28 


7 


14:9 


47'41" 


2.3 


0.4 


17 


8 


14*7 


47'44" 


1.3 


0.5 


10 


9 


1X6 


47'48" 


1.8 


0.3 


18 


10 


1415 


47'52" 


0.8 


0.3 


160 


11 


15!3 


47'51" 


2.2 


0.3 


8 


12 


15*3 


47'58" 


2.1 


0.3 


13 


13 


14:6 


48'05" 


0.6 


0.3 




14 


15!3 


48'06" 


1.8 


0.3 


10 


15 


1K6 


4813" 


2.1 


0.3 


178 


16 


13?5 


48'15" 


1.8 


0.6 


19 
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Table 5. Reflection Nebulosity in W42 



Patch 


RA J2000 Decl. 


Pk{%) 


M°) 




18 fc 38 m 


-6° 






1 


14!69 


47'54" 


13 


100 


2 


14*66 


47'57" 


32 


97 


3 


14!55 


48'02" 


20 


97 


4 


14!50 


48'06" 


20 


106 


5 


14*47 


48' 11" 


26 


101 



Table 6. Stars in Cepheus A 



Star RAJ2000 Decl. H — K P H (%) e P (%) 6 H {°) 
22 fe 56 m +61° 



1 


16*0 


02'19" 


2.75 


5.9 


0.2 


123 


2 


1753 


02'24" 


0.73 


0.9 


0.2 


154 


3 


19!8 


02'21" 


0.61 


0.8 


0.2 


54 


4 


22!5 


02'01" 


1.39 


4.9 


0.3 


63 


5 


22!6 


01'42" 


1.45 


2.1 


0.3 


94 
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